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Abstract  This  paper  discusses  on  the  advantages  and 
shortfalls  of  the  integration  of  the  2  navigation  systems,  namely 
the  GPS  and  the  INS,  and  then  it  discusses  a  method  of  integrating 
them  by  use  of  the  pseudo  range  and  pseudo  range-rate  measurements 
of  the  GPS-receiver  to  integrate  GPS  with  INS.  Because  the 
measurement  errors  of  the  GPS  are  all  different,  the  numerical 
simulation  of  the  integrated  navigation  capability  of  a  plane  had 
to  be  computed  under  different  conditions  of  the  differential 
machines,  and  details  of  the  accuracy  analyses  were  also  made 
accordingly.  The  results  showed  that  such  integrated  system  of 
navigation  provided  the  accuracy  which  could  not  be  achieved  by 
either  the  GPS  or  the  INS  alone.  By  integrating  the  differential 
GPS  receiver  with  the  INS,  the  errors  in  position  could  be  further 
reduced.  The  integrated  navigation  thus  opens  up  a  new  area  of 
application  to  be  used  for  a  plane  to  approach  an  airport  for  a 
landing. 
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I.  INTRODUCTION 


The  inertial  navigation  system  of  all  the  aircraft 
facilities  has  a  unique  function  to  provide  the  accurate 
communication  information.  However#  the  errors  in  position 
will  increase  as  time  goes  by.  To  overcome  such  shortcoming# 
often  other  navigation  facilities  installed  in  the  aircraft 
are  being  used  in  tandem  in  order  to  obtain  informations  of 
most  superior  quality#  statistically  speaking.  The  GPS  as 
the  global  positioning  system  can  make  accurate 
determinations  of  3-dimensional  positions#  speed  and  time# 
but  it  cannot  provide  the  needed  information  of  the  plane# 
such  as  the  signal  to  show  the  altitude  of  the  plane  or  the 
effects  of  the  plane  movements#  in  order  to  have  its  motion 
tracked  or  traced.  By  integrating  the  INS  with  the  GPS#  one 
can  overcome  this  shortcoming#  and  form  an  ideal  integrated 
system.  Furthermore#  other  traditional  functions  of  the  GPS 
remain  intact  to  continue  its  reconnaisance  work  and  by  use 
of  the  differential  method  or  capability  of  the  ground 
surface  pseudo  position  to  improve  the  reliability  on  the 
accuracy  of  the  position  determination  by  the  GPS. 
Consequently,  analyzing  the  accuracy  of  the  integrated 
navigation  system#  making  accurate  estimations  and  utilizing 
the  integrated  system  for  some  new  functions  all  have 
important  significances . 
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II.  THE  INTEGRATION  METHOD  OF  GPS/ INS 

The  integrating  method  basically  can  be  classified  into 
2  kinds  : 

1.  By  use  of  the  GPS  receiver  (below,  to  be  just  called  as 
the  GPS  instruments)  to  obtain  the  information  about 
the  position  of  a  plane  and  its  measured  speed,  and  to 
integrate  them  with  INS  (1). 

The  information  sent  out  by  the  plane  can  be  the 
preliminary  data  to  compute  its  position  and  speed  with 
respect  to  4  stars  (such  method  is  commonly  known  as  the 
point-solution  method),  or  it  can  be  the  estimated  values  of 
the  position  and  speed  by  use  of  "Carmey"  filter  wave 
machine.  The  advantage  of  integration  is  to  provide 
additional  navigational  information,  to  simplify  the 
measuring  procedure,  and  to  reduce  the  filter  wave 
computation.  But  because  the  measuring  equations  are  unable 
to  describe  accurately  the  statistical  and  the  correlative 
characteristics  of  the  noises  relevant  to  the  position  and 
velocity,  and  because  of  the  effects  on  the  filter  wave,  the 
influence  becomes  even  more  prominent  due  to  the  difference 
of  the  conditions  in  the  satellite  geometrical  distribution. 

2 .  Utilization  of  the  Pseudo  Range  aDd  Pseudo  Range-Rate 
of  the  GPS  instrument  to  integrate  with  the  INS 

The  pseudo  range  and  pseudo  range-rate  could  be  the 
actually  measured  values  (1,2)  or  they  could  be  the 
estimated  values  obtained  from  "Carmey"  filter  wave  machine 
(3).  From  the  measuring  equations,  one  can  directly 
describe  the  errors  of  the  pseudo  range  and  pseudo 
range-rate  to  improve  the  integration  and  the  result  is  less 
affected  by  the  difference  in  the  satellite  geometrical 
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distribution.  But  the  parameter  computation  in  the 
measuring  equations  is  complicated  and  the  computed  values 
are  relatively  large. 

III.  SIMULATIONS  OF  THE  NAVIGATIONAL  FUNCTIONS  OF  THE 

INTEGRATION  OF  THE  GPS  WITH  THE  INS  BY  UTILIZING  THE 
PSEUDO  RANGE  AND  PSEUDO  RANGE-RATE 

By  use  of  pseudo  range  and  pseudo  range  rate,  an 
example  of  the  integration  was  worked  out  to  simulate  its 
functions  ;  the  object  to  simulate  is  the  integration  of  the 
plane  platform  type  of  the  INS  of  an  average  accuracy  which 
goes  through  the  ground  as  the  initial  step  and  the  GPS 
instrument  of  a  single  channel-multiple  route  reusable  C/A. 
The  INS  navigation  system  assumes  the  horizontal  coordinate 
system  (the  W-system)  with  a  floating  azimuthal  direction 
for  its  coordinate  system ,  the  geographical  coordinate 
system  (the  g-system)  as  the  error  equation  coordinate 
system  (4),  and  the  star  base  of  the  GPS  places  the  entire 
distribution  of  18  stars  under  its  consideration. 

1.  Introduction  of  simulation 

(1)  Integrating  the  states  of  filter  wave.  There 
are  altogether  28  states  which  include  18  error-states  of 
the  INS  of  the  position  errors  ( SL ,  SX,'  Sh)  or  (SR.,  SRwrSh). 

velocity  errors  ( SVH ,  5VW,  5VZ),  ,  angular  errors  of  the 

plane  platform  (<&„,  Qw,  bias  in  the  acceleratometer 

(  ix'  '  gyroscopic  bias  (thX,  tkr,  t^)  and  the 

gyroscopic  floating  displacements  in  one  Markovian  Process 
(emX,  tmY,  emZ);  10  error-states  of  the  GPS  instrument  which 

include  6  clock  errors,  phase  error  (SQ).  frequency  error 
iSJl,  frequency  variation  error  (Srj) , ;  and  frequency  error  due 
to  the  velocity  of  the  plane  (£f;);  ;  4  bias  in  pseudo  range 

measurements  ( Spbi ,  /=1,..,4).  .  Among  them  the  average 
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values  of  «• ,  e r*  and  (i**X,Y,Z,  which  are  the  axial 

directions  of  the  W-system)  turned  out  to  be  zero  and  the 
average  square-root  errors  of  them  are,  respectively, 


'  (£{  *})*  -  0.5  x  10  *g,  (£{efj)i  =0.0!(  °  )/  h, 
jf£{0>'  =  0.005(  °  )/  h  f  ... 


r  i 


For  the  INS  error  state  equation  and  the  GPS  instrument 
clock  error  equation,  see  references  (4,  5). 

(2)  The  Measuring  Equations.  If  the  reference 
system  is  the  g-system,  the  correlative  formulas  of  pseudo 
range  and  pseudo  range-rate  are,  respectively, 

",  -K*.  -•*■„>’ +(r.  ~  +s,v  I 

",  +  -  rjY'  -  r„)  +  (z,  -z„xz, -z„>] 

+  <54>  +  6paj 

where  (XV  .  Z  g.)  and  (x f.,  Ff  .  ^.f)'  are  the  coordinates  and 

the  rate  changes  of  the  coordinates  of  the  carrier  body  with 
respect  to  the  g-system,  (Tj,  Y^,  Zj'  and  '  (x^,  y^, 
are  the  coordinates  and  the  rate  changes  of  the  coordinates 
of  the  j-th  star  with  respect  to  the  g-system,  ‘/ry  is  the 
relative  distance  between  the  j-th  star  and  the  carrier,  and 
5pwj  and  SpwJ  are,  respectively,  the  white  noises  of  py |  and 
p  By  use  of  the  INS,  the  navigation  coefficents  Le, \  1,. 

hc  and  rw<.  V  Wt„  V Ze  can  be  computed,  and  from  the 
trajectory  data  of  the  satellites  in  the  telecommunication 
messages  of  the  GPS,  one  can  obtain  the  computed  distance 
PM,  .  If  the  measured  value  of  the  j-th  star  is 


0) 

(2) 


Z j  not  only  has  bearings  with  the  GPS  errors  , 

5PWJ>6p9J)'  /  but  it  is  also  related  to  the  INS  errors  { SL , 

nSXMJSV fltSVw,SV2Y  ,  and  thus  from  it  one  can  find  out 

the  measuring  equation  which  nas  bearings  with  the  state  of 
the  filter  wave  device. 

If  the  frequency  of  the  simulated  filter  wave  is  2s, 
following  the  procedure-treating  computations  and  by  use  of 
the  method  shown  in  reference  (6),  one  can  choose  4  stars  as 
the  useful  ones. 

(3)  The  orbit  of  the  plane  chosen  by  Simulation. 
The  orbit  of  the  plane  is  in  5  stages.  The  plane  was 
accelerated  for  200  s  (a  =  1.5  m/s1))/  it  climbed  up  for  100 
s,  flew  horizontally  at  speed  of  300  m/s  for  3,400  8,  then 
making  a  180°  turn-around,  it  flew  horizontally  back  home. 
The  total  flight  was  about  for  7,500  s. 

2.  The  Results  of  Simulation 

By  use  of  pseudo  range  to  compute  the  noise  and  bias, 
the  following  3  kinds  of  integration  simulation  calculation 
were  made  : 

(1)  Simulation  1.  In  simulating  the  basic  error 
data,  namely  the  average  square-root  errors  of  S^w,,  5^  and 
5,.'  ,  their  values  are  as  follows  (1,7,8): 

<5 .  =  0.05m /s,  <5  =  10m.  6  =  10m.  \ 

fw  f*  rw 
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Fig 


.  1  Simulation  1  The  propagation  curve  of  velocity 

errors 


Fig.  2  Simulation  1  The  propagation  curve  of  shape  errors 


The  average  square-root  deviation  of  velocity  errors , 
attitude  errors  and  position  errors  from  the  results  of 
simulation  (namely  the  average  square-root  deviation  of  the 
errors  after  the  corrections)  as  well  as  their  propagation 
curves  are  all  shown  separately  in  Fig.l  -  Fig. 3  ;  the 
geometrical  accuracy  factor  GDOP  of  satellites  and  the 
expansion  coefficients  of  the  position  errors  of  Northwest 
Sky  NDOP/  WDOP  and  VDOP  propagation  curves  are  shown  in  Fig. 
4.  The  abrupt  change  in  Fig.  4  was  due  to  the  replacement 
of  satellites. 


(2)  Simulation  2.  When  selective  availability 

(SA)  was  assumed  for  the  C/A  ratio  in  America,  the  bias  in 
pseudo  range  increased.  Simulation  2  chose  =  30  m  (8), 

and  simulation  was  affected  by  SA.  The  results  of  the 
simulation  were  similar  to  those  from  simulation  1,  if 
position  errors  were  increased  slightly. 

(3)  Simulation  3.  A  simple  recognition  of  the 
differential  GPS  could  entirely  compensate  for  bias  of  the 
pseudo  range.  Therefore  if  =  0,  the  integration  between 
the  simulated  differential  GPS  and  the  INS  was  effective. 

As  compared  to  Simulation  1,  the  position  errors  went  down 
8 lightly. 
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IV.  THE  ACCURACY  ANALYSES  OF  THE  SIMULATED  RESULTS 


1.  From  the  results  of  simulation  1  and  simulation  2  to 
estimate  the  accuracy  analyses  of  position  errors  and 
velocity  errors,  the  estimated  accuracy  of  position  errors 
and  velocity  errors  turned  out  to  be  pretty  high,  and  the 
reasons  seem  to  be  as  follows  : 

(1)  The  results  of  a  direct  integration  of  the 
highly  accurate  GPS  pseudo  range  and  pseudo  range  rate  with 
either  the  position  error  or  the  velocity  error  of  the  INS. 
By  use  of  a  simplified  situation  one  can  analyze  the 
declining  rate  of  the  INS's  position  errors.  Suppose  a 
system  is  integrated  only  with  the  pseudo  range,  and  from 
the  GPS  instrument  one  can  directly  measure  the  3 
coordinates  of  the  distance  between  a  satellite  and  the 
carrier  in  the  g-system.  The  measurement  error  is  the  white 
noise,  and  the  average  square  root  error  is  (i  *  N,  W, 

Z).  If  one  neglects  the  clock  error  of  the  GPS  instrument, 
and  under  the  condition  of  non-diagonal  elements  of  the 
estimated  average  square  root  error  matrix  P  (k-1)  of  the 
filter  wave  device,  one  can  obtain  the  approximate  formula 
for  the  estimated  average  square  root  error  PIRi(k) 

(i=N,W,Z)  of  the  position  error  at  time  k  as  follows  : 

r,  /(/>„,»/*  - !)+«’)  (4) 

The  above  expression  shows  that  •  Since  the 

frequency  of  the  filter  wave  in  general  can  be  assumed  to  be 
1-58,  and  because  P m  (k/k-1)  cannot  be  too  much  greater 
than  Pm  (k-1),  one  can  approximately  recognize 

.  •>:  Pm  (*/*-!)  =  <>/»,„  .(*-!> 
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Indeed  Ptu  (k)  •  After  several  estimations, 

one  can  let  PtK,  (k)  *  The  reason  is  this.  If 

the  estimated  state  is  in  the  condition  capable  of  being 
directly  measured,  "Carmey"  filter  wave  utilizes  all  the 
timely  measured  values  and  thus  it  can  filter  out  all  the 
noise.  Also  it  is  because  as  shown  in  Reference  (7)  the 
filter  wave  equipment  can  filter  out  2/3  of  the  measurement 
noise. 

(2)  The  results  of  double  integration  of  the  pseudo 
range  and  pseudo  range  rate  of  the  GPS  with  both  the 
position  error  as  well  as  the  velocity  error  of  the  INS. 

By  use  of  the  GPS  pseudo  range  to  integrate  with  the  INS,  it 
not  only  improved  the  position  error  but  because  the  pseudo 
range  also  has  the  capacity  to  observe  the  velocity  error  it 
could  improve  the  velocity  error  as  well.  By  use  of  the  GPS 
pseudo  range  rate  to  integrate  with  the  INS,  beside 
improving  the  velocity  error,  it  could  also  suppress  the 
increase  of  the  position  errors.  To  prove  the  above 
mentioned  analyses,  a  supplementary  simulation  was  carried 
out  on  the  integration  only  by  the  pseudo  range.  The 
results  showed  that  accuracy  in  the  estimation  of  the 
position  error  had  crone  down,  but  for  the  velocity  error  it 
definitely  had  the  effectiveness  of  the  estimation,  and  thus 
it  explains  that  the  double  integration  can  lower  both  kinds 
of  errors  further  down  by  one  step. 

<3)  The  integration  enables  the  pseudo  range  bias 
to  get  estimated.  In  analyzing  the  observability  of  the 
system,  one  can  see  that  the  pseudo  range  bias  is  an 
observable  quantity.  The  simulated  results  showed  that  as 
the  times  of  estimations  increased,  the  average  square  root 
of  the  pseudo  range  bias  went  down  from  10  m  to  around  3  m, 
while  bias  of  30  m  could  go  down  to  7  -  8  m.  This  is  a 
great  effect  in  suppresing  the  measurement  errors  of  the 
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pseudo  range  extensively.  Such  quality  to  rely  heavily  on 
the  GPS  instrument  is  not  appropriate. 

(4)  Integration  can  suppress  unfavorable  effects 
of  the  satellite  geometrical  distribution.  As  widely 
known,  one  can  use  GDOP  to  describe  the  influence  of  the 
geometrical  distribution  of  satellites  on  the  point-solution 
method  in  making  position  determinations.  Below,  analysis 
will  be  made  with  one  simplified  situation.  Suppose  that  a 
satellite  is  on  the  same  horizontal  plane  as  the  carrier, 
and  one  can  neglect  the  clock  error,  just  assuming  the  GPS 
measurement  error  as  a  white  noise  and  the  average  square 
root  deviation  .  Now  the  pseudo  range  of  the  GPS  is 
integrated  with  the  INS  and  then  the  simplified  equation  for 
the  position  error  of  the  INS 

s*.-*.  (5> 

where  W H  and  W w  are  the  white  noise.  The  distances  of  the 
carrier  to  the  2  stars  were  ( p,  and  p, )  while  those  computed 
by  the  INS  were  (p(|  and  pfJ),  and  the  measured  values  were 
the  differences  between  these  2  sets  of  values;  that  is. 


where  6J  (J-U  2)  is  the  azimuthal  angle  of  the  line  which 
connects  the  j-th  star  to  the  carrier.  To  approximate  ♦»/ 
(I  is  the  unit  matrix),  it  becomes  the  transformation  matrix 
after  dispersing.  Neglecting  estimated  average  square  root 
error  matrix  P(k-l)  non-diagonal  elements  and  going  through 
the  solution  of  the  Carmey  filter  wave  equation,  one  can  get 


the  average  square  deviation  of  the  horizontal  errors  at 
time  k  as  follows  : 


/>*W 


[  1  +  f Is 


COJ0  •" 


where  P1  (*X«  =*  1,2)  is  the  i-th  diagonal  elements  of  P(k) 

matrix.  And  with  0**  (5, -QJ  t-90*. 

,  ..  r  ]*  , |  j.  |  i _ 

'»  W,*,  Ll+Cf.  +  ^J  '  ” 

[P((l  +ji>i10  —  sin*8l  —  ««,S1)  + [Pj(l  —  coti0i  ~cos,8J)]  (8) 

where  J®,  and  P,  are  2  diagonal  elements  of  P  (k/k-l)  In 
equation  (7),  the  factor  VJ / cos$  is  the  geometrical 
accuracy  factor  of  the  2-dimensional  position  determining 
system.  When  0,-0,*=  90°  ,  namely  0  =  0,  the  geometrical 

distribution  is  the  best  ;  that  is# 

V«-o)-/>ll+^ (9) 

*(0-0) -I" 1  V  mK^  <  1  (10) 

rir! 

When  one  looks  at  the  contents  of  R,  one  sees  that  it 
represents  “Carmey"  filter  wave  which  was  at  the  estimated 
base,  namely  the  k-th  estimated  effect.  When  0=0,  the 
estimated  effects  and  the  2  horizontal  directions  (N,  W)  all 
had  the  same  estimation  effects  from  the  independently 
measured  values.  When  0=£O,  if  the  approximation  allows 
Pi^Pt,  ,  then  from  formula  (8)  and  (9),  Pt>Pjn in,  and 
always  .  To  explain  ,  "Carmey"  filter  wave  can 

also  suppress  any  unfavorable  influence  from  the  satellite 
distribution.  A  similar  analysis  on  the  integration  with  the 
velocity  reached  the  same  conclusion. 


Although  the  above  analysis  was  made  by  a  simplified 
2-dimensional  situation,  yet  the  results  can  be  extrapolated 
to  the  actual  integrated  situations  of  the  GPS  and  the  INS. 
The  result  of  simulations  confirmed  such  characteristics; 
that  is,  although  the  velocity  error  and  position  error  both 
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receive  the  GDOP  influence,  but  the  influence  is  relatively 
small.  Below,  the  accuracy  of  the  GPS  instrument  navigation 
and  these  of  various  integrated  navigations  are  all  shown  in 
the  table  for  comparisons.  Among  them,  there  are  2  kinds, 
namely  the  GPS  instrument  navigation  point- solution  method 
and  the  filter  wave  solution  method.  In  the  point-solution 
method,  the  coefficient  of  expansion  HDOP  and  VDOP  was  used 
to  calculate  the  navigation  accuracy,  while  the  filter  wave 
solution  method  filtered  out  2/3  of  the  measurement  noises, 
and  by  utilizing  this  effect  it  can  calculate  the  navigation 
accuracy.  The  integrated  accuracy  used  the  last  stage  of 
the  simulated  flight  as  its  accuracy.  It  also  used  HDOP  - 
1.6,  VDOP  =*  2.7,  re$  -  10  m  (30  m),  <Vv  =  10  and  33 

0.05  m/s.  From  the  table,  one  can  easily  see  the  effects  of 
improvement  in  the  navigation  accuracy  due  to  integrations. 


Table 


(1)  Classifications  of  Errors _ position  Errors  (m) 

Velocity  Errors 

(2)  Different  kinds  of  Integrations  and  Solution  Methods 
Average  Square  Root  Deviation  Horizontal  Vertical 
Total  Error  Horizontal  Vertical  Total  Error 

(3)  The  GPS  instrument  point-solution  method 

The  GPS  instrument  filter  wave  solution  method 


The  differential  GPS  instrument  filter  wave  solution  method 


Integration  of  Basic  GPS/INS  (Simulation  1) 

Integration  of  GPS  with  large  bias/INS  (Simulation  2) 
Integration  of  differential  GPS/INS  (Simulation  3) 

2.  Attitude  Error  Estimation  Accuracy.  When  the  INS 
is  integrated  with  the  GPS  at  the  initial  condition,  it  can 
only  suppress  the  increase  of  the  horizontal  error-angle 
with  respect  to  the  plane  platform.  However,  there  always 
exists  the  space-to-air  aiming  capacity  (9)  in  such  kind  of 
integrations,  and  consequently  the  estimating  action  of  this 
kind  of  integration  under  the  flight  condition  of  an 
airplane  toward  the  azimuthal  error  is  clear.  To  explain 
Fig.  2,  going  through  the  horizontal  acceleration  and  the 
180°  Turn-around,  the  azimuthal  error  would  decline  from  3' 
to  0.5  -  0.7' ,  which  is  very  good. 

3.  The  integration  unequally  improves  the  original 
errors  of  the  INS  and  the  clock  errors  of  the  GPS 
instrument.  In  them,  using  the  vertical  acceleration  to 
compute  bias  and  making  an  estimation  of  the  error  which  is 
sensitive  to  the  clock  acceleration  are  especially 
noticable.  The  reason  is  obvious.  Improving  the  original 
errors  is  more  advantageous  when  the  GPS  signal  is 
interrupted,  and  in  this  case  the  functions  of  the  INS  are 
improved . 

It  needs  some  explanation  if  one  considers  the  error  in 
the  movement  of  the  GPS  instrument  or  when  one  uses  a 
stepping-down  wave  machine.  The  above-mentioned  filter  wave 
effect  receives  a  constant  influence. 
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V. 


CONCLUSION 


1.  When  the  GPS  instrument  made  a  double  integration 
with  both  the  position  and  velocity  of  the  INS,  such 
integration  made  the  bias  of  the  GPS  pseudo  range 
observable,  and  it  could  also  suppress  the  effects  of  GDOP, 
and  therefore  the  estimation  effect  utilizing  the 
integration  of  the  pseudo  range  and  pseudo  range  rate  of  the 
GPS  with  the  INS  was  good. 

2.  The  integration  allowed  the  INS  to  possess  the 
space-to-air  aiming  capability.  Thus  when  a  plane  was 
accelerated  and  turned  around,  it  could  clearly  make  an 
estimation  of  the  INS  azimuthal  error-angle.  The  results 
from  simulation  explained  that  the  error  could  go  down  to 
below  1 ' . 


3.  After  the  GPS  instrument  integrated  with  the  INS 
of  only  an  average  accuracy,  the  position  error  of  the  last 
stage  of  the  flight  according  to  the  simulation  was  already 
about  7  m.  If  the  differential  GPS  was  used  to  integrate 
with  the  INS,  the  error  could  even  go  down  to  around  3  m, 
and  moreover  the  main  component  of  the  errors  was  the 
vertical  error  component  and  therefore  if  the  system 
expressed  altitude  in  the  integrated  radio,  the  total 
position  error  could  even  be  reduced  at  lower  altitudes. 
Such  kind  of  position  accuracy  is  what  the  aircraft 
navigation  craves  for.  Consequently,  this  integration 
system  can  be  expected  to  open  up  a  new  application  area  to 
let  airplanes  make  automatic  entry  to  an  air-field  for  a 
landing. 
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